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Abstract — A diffusion-based molecular communication system 
has two major components: the diffusion in the medium, and the 
ligand-reception. Information bits, encoded in the time variations 
of the concentration of molecules, are conveyed to the receiver 
front through the molecular diffusion in the medium. The 
receiver, in turn, measures the concentration of the molecules 
in its vicinity in order to retrieve the information. This is done 
via ligand-reception process. In this paper, we develop models 
to study the constraints imposed by the concentration sensing 
at the receiver side and derive the maximum rate by which a 
ligand-receiver can receive information. Therefore, the overall 
capacity of the diffusion channel with the ligand receptors can 
be obtained by combining the results presented in this paper 
with our previous work on the achievable information rate of 
molecular communication over the diffusion channel. 

I. Introduction 

Molecular diffusion plays a key role in different forms 
of communication in nature. At the microorganism scale, 
molecular signals are used for communication and control 
among cells. For example, one of the well-known communi- 
cation primitives among cells is a phenomenon called Quorum 
Sensing. Quorum sensing is a communication process among 
cells in which bacteria monitor the population density of their 
type (and possibly other types of bacteria) via the production 
and detection of special molecules. Quorum sensing allows 
bacteria to synchronize the behavior of the group, and thus 
act as a unit fTll-l3ll. 

The communication can also appear in other forms such as 
pheromone communication and calcium signaling. Pheromone 
communication is performed using specific type of molecules 
released by plants, insects and other animals which trigger 
specific behaviors among the receptor members of the same 
species. Pheromones of a particular type, after being released 
into the environment, diffuse in the air until they are captured 
by the same receptor type. Cells absorb calcium molecules in 
response to various stimuli that open/close particular channels 
on the cell membrane. The molecular information in the 
variation of the calcium ions concentration is propagated 
both inside and outside the cell, causing a variation in the 
electrical charge of the cell membrane and, subsequently, the 
transduction of the information into an electrical signal. 

Due to various limitations such as size, arguably the dom- 
inant form of communication in nano scales is via molecular 
signaling, which is fundamentally different from conventional 
electromagnetic-based communication. Recent advances in 
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bio-nano technology has motivated research on designing 
nanoscale devices to perform tasks similar to their biological 
counterparts. This opens up several potential applications. 
For example, it could lead to molecular based networks that 
can be deployed over or inside the human body to monitor 
glucose, sodium, and cholesterol, to detect the presence of 
different infectious agents, or to identify specific types of 
cancer. Such networks would also enable new smart drug 
administrative systems to release specific drugs inside the body 
with great accuracy and in a timely manner. To enable all such 
applications, the communication is the key. 

The wide range of potential applications has inspired de- 
velopment of new theoretical frameworks for molecular com- 
munication 0-0. The focus of these studies is to encode 
the information into the time of the release of molecules 
while the propagation from the transmitter to the receiver is 
governed by Brownian motion (without drift [5| and with 
drift El, ||6l). Although Brownian motion of particles is 
basically the same process as diffusion, natural organisms 
do not look directly at the Brownian motion of a large 
particle. Instead, they measure diffusion of small molecules by 
following changes in concentrations. Therefore, we consider 
an alternative molecular communication system governed by 
diffusion process in the medium, referred as the diffusion 
channel (8). Arguably, the most dominating form of the 
communication at the micrometer scale is diffusion based 
molecular communication, i.e., embedding the information in 
the alteration of the concentration of the molecules and rely 
on diffusion to transfer the information to the destination. 
The communication among bacteria, calcium signaling, and 
many others can be reduced to diffusion based molecular 
communication. Another body of work in this field involves 
the study of the Quorum Sensing as a network and mapping 
the Quorum Sensing to consensus problem under diffusion- 
based molecular communication 0. 

In our previous work JS], we studied the simpler and more 
practical discrete systems. There, we studied the achievable 
rate over the diffusion channel when we are limited to the 
amplitude modulation (on-off keying), i.e., a discrete diffusion 
channel. One of the important problems that where left open 
in (H was the reception process. As explained later, the 
dynamics of the reception in molecular communication is 
fundamentally different from the conventional receivers. As 
a result, the receiver can limit the capacity of the molecular 
communication. In this work, we investigate the dynamics of 
the molecular receptor and its effect on the capacity of the 
molecular diffusion communication. Further, using the results 



in this paper, we verify some of the basic assumptions we 
made in our previous work. 

We consider a communication scenario consisted of a 
transmitter and a receiver communicating via the diffusion 
channel, as shown in Fig. Q] In this context, the transmitter is 
capable of releasing molecules to the medium to change their 
concentration according to information bits [8|. Similarly, the 
receiver is capable of absorbing the molecules or chemical 
signals, for example, by using ligand-receptor binding which 
is a trans-membrane receptor protein on a receiving cell [ITOll . 
fTT1l . Ligand-receptors provide a bin for molecules to bind to. 
This process triggers a signaling inside the cell indicating the 
reception of the molecule. The receiver has a large number of 
binding places using which it can estimate the concentration 
by averaging over all binding places. The performance of the 
ligand-receptor has a profound impact on the reception of 
information. Once a molecule is bound to a receptor, it takes 
some time for the receptor to reset, we later address this in 
more detail. Our goal is to find the capacity in such a diffusion- 
based molecular communication channel using the models we 
develop for the ligand-receptor. 

The rest of the paper is organized as follows: In Sec. HH 
we review the model for diffusion-based molecular commu- 
nication. In Sec. HID the model for an ideal receiver and its 
capacity is studied and then it is extended to a Markov-based 
model for receptors. Sec. |IV] addresses the link between this 
paper and our previous work and concludes the paper. 

II. Molecular Communication Model 

The paradigm for diffusion-based molecular communication 
is depicted in Fig. [T] We assume a source with binary 
alphabets. The information bits are encoded into alteration 
of the rate of molecule production by the transmitter T. 
In particular, we consider pulse-shaped modulation, i.e., the 
transmitter either releases molecules over a time period or 
remains silent. The emitted molecules are transferred through 
the medium by diffusion process and reach the receiver R. The 
flow of molecules from the regions with high concentration 
(around T) to the regions with lower concentration (around 
R) is governed by the diffusion equation. The second Fick's 
law of diffusion, as discussed in E), gives the concentration of 
molecules at each position in the medium over time, proportion 
to the input rate. We can view the diffusion in the channel 
as a low-pass linear filter which takes the production rate as 
the input and gives the concentration at the receiver as the 
output via convolution with the impulse response (8). The 
low-pass frequency response of the diffusion channel restricts 
the rate of the information exchange; the receiver does not 
observe the high frequency elements which in turn means 
that the rapid changes in concentration of molecules at te 
receiver are absent. In JS], we also observed that there is an 
inherent memory in the medium. Molecules tend to linger in 
the medium after each transmission and some amount of time 
is needed for the molecules to diffuse away and channel to 
be reset. In (8), we considered all the above constraints and 
obtained the capacity of the diffusion channel; assuming a 
pulse amplitude modulation (of the molecule production). 
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Fig. 1. Basic molecular communication model: the transmitter, diffusion 
medium, and the receiver. 

Once the diffusion process transfers the molecules to the re- 
ceiver, the receiver R measures the concentration of molecules, 
using a set of binding sites that become active once they trap 
a molecule. Each binding site can trap one molecule at a time 
and cannot accept any more until the reaction to the previously 
trapped molecule is completed. The model for the receptor 
is shown in Fig. [2] As the concentration p increases, the 
probability p that a molecules is trapped and hence the number 
of active sites increases. This way, the receiver can estimate the 
concentration in the medium. Note that the receiver samples 
these binding sites (for being in active state) synchronously at 
the end of each bit transmission (as pulse-shape keying PSK 
in the conventional communication systems). The number of 
active binding sites will determine the received concentration. 
The information can then be extracted from the variation of 
the concentration of molecules at R. 

As mentioned earlier, the entire diffusion-based molecular 
communication system has two major components: the diffu- 
sion in the medium, and the ligand-reception. The capacity 
due to the diffusion channel is obtained in [8]. In this paper, 
we study the fundamental limits on the information rate due to 
the ligand-reception. The overall capacity of the system, which 
is the maximum conveyed information per channel use, can 
be easily obtained by having the capacities of the individual 
components. 

III. Modeling the Ligand Receptors 

Assume a receiver equipped with N binding sites/receptors. 
In every moment, a receptor is either active (a molecule 
is trapped in it) or inactive (empty). We assume all the 
receptors are identical, observe the same concentration, and act 
independently. The process of determining the concentration of 
molecules at the receiver is as follows: the molecules arriving 
from the channel make up a time varying concentration p(t) 
at the receiver vicinity. The higher the concentration around 
the receiver, the higher the probability that a molecule gets 
trapped in a receptor. We denote by p the probability that a 
single molecule being bound to an empty ligand-receptor when 
the concentration p(t) = p. In other words, 

p = P (Binding | p & the receptor is empty) (1) 
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Therefore, the problem is to obtain fp{p) that achieves the 
capacity in ©. 

The capacity of the receiver depends on the model assumed 
for functionality of the receptors. In the following, we study 
two models to describe the process of estimating the con- 
centration of molecules. First, the ideal receptor model is 
discussed. Then, we extend this model to capture the more 
complicated functionality of the receptors using a Markov 
chain model. 



Fig. 2. Illustration of the ligand receptor mechanism for sensing the 
concentration. 



The function that maps the concentration of molecules to the 
probability of the absorption is assumed to be a monotonically 
increasing one, [0, oo] —> [0, 1]. As discussed later, p G [0, 1] 
is a realization of a random variable P. 

To each receptor, we associate an indicator random variable 
Xi, i G {1,2, ...,7V}, where X^s are Bernoulli distributed 
with parameter p, i.e., P(Xi = 1) = p. Let 1^' be the vector 
containing the outputs of all the receptors. The vector X^ is 
the input used to estimate p. It is known that Y^iLi Xi is the 
sufficient statistics for estimation of p. This implies that the 
receiver functionality is to add up the values of all receptors at 
every sampling instant. Let P be the best unbiased estimator 
for p. Since P = YliLi Xi, the expected value and variance 
of P, given p, can be easily derived as E[P] = p and Var[P] = 

&>(i-p)- a 

Clearly, P depends on the value that P takes. In other 
words, we have 
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We use the shorthand P(i\p) for P (p = jj\P = p) in ©. 

Hence, by computing P, the receiver has an estimate for p 
which in turn is a function of p. The receiver maps back the 
estimate P to a level of concentration p. Thus, we have a 
chain of processes as shown in Fig. [2] We observe that except 
for the relation between P and P, all the other relations are 
deterministic and one-to-one. Therefore, to study the capacity 
of the ligand receptor, it is sufficient to consider only the 
mutual information between P and P. This implies that 
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where C is the capacity due to the ligand reception and fp (p) 
is the distribution function of P. In addition, H(P) is the 
entropy of discrete random variable P which takes values from 
the set {-^, i G 0, 1, iV}. Using ® and the definition of 



A. Ideal Receptors 

An ideal receptor is defined as a receptor which is always 
ready to receive molecules. In this model, molecules are 
assumed to occupy the receptors for a negligible amount of 
time. Hence, we assume that the molecules leave the receptor 
instantly after absorption and the condition on p that needs the 
receptors to be empty is always satisfied. This means that the 
parameter P would be directly an estimate for the probability 
in (Q]). Therefore, 



(P;P)=/(P;£ti^ 
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In this case, the problem of finding the capacity achieving 
distribution in © resembles the corresponding problem in 
source coding which is discussed in |[T2ll . There, the problem is 
to find the distribution that maximizes the mutual information 
between an observed sequence and the unknown source pa- 
rameter 0. It is shown that the capacity achieving distribution 
is Jeffreys Prior which is proportional to the square root of 
the determinant of the Fisher information X: 
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In our setup, we intend to maximize the mutual information 
between the sequence 1^' and the unknown parameter p. 
By comparison to the source coding problem we have: 

Lemma 1. The capacity achieving distribution on P that 
maximizes the mutual information I(P,X^) follows the 
Jeffery's prior [12]. 

It is easy to verify that the Fisher information contained in 
N independent Bernoulli trials with parameter is 

N 



1(0) 



1 



Therefore, the Jeffrey's prior for a Bernoulli random variable 
is given by the Arcsine distribution: 
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This capacity achieving distribution is shown in Fig. [3] The 
resulting distribution for P shows that, in order to achieve 
the capacity, the transmitter should produce molecules in a 
way that approximately, makes P close to either 1 or 0, in 
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Fig. 3. Capacity-achieving distribution for p 
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Fig. 4. Capacity versus the number of receptors 



each transmission. This would also justify the use of a binary 
level pulse transmitter in our previous work [0. Intuitively, 
when p is chosen to become closer to ~, the variance p(l—p) 
of the Bernoulli output of each receptor is increased. This 
would reduce the mutual information in ©. The numerical 
results for the capacity of the ligand receptors versus the 
number of receptors N is shown in the outermost plot in 
Fig. HI Expectedly, the capacity is monotonically increasing 
with respect to N. 

B. Markov Chain Model for Receptors 

In fT3ll . a more realistic model to address the lingering of 
molecules in receptors is used. In order to improve our ideal 
receptor model, we consider a Markov Chain (MC) model to 
better capture performance of receptors. This is due to the fact 
that a molecule which is trapped by a receptor does not leave 
instantly and will occupy the binding site for a random amount 
of time. During this interval, the receptor remains occupied 
and no other molecule is able to bind to the receptor. Hence, 
the probability that each receptor accepts a molecule is no 
longer p. Instead, some randomness would exist because of 
the non-ideal performance of the receptors. The MC model 
for each receptor is shown in Fig. [5] In this figure, state 



Fig. 5. Markov chain model of a ligand receptor 

corresponds to an empty receptor, i.e., X{ = 0, whereas state 
1 corresponds to the receptor with a trapped molecule, i.e., 
Xi = 1. Further, p is the same as the previous section defined 
in [TJ and q is the probability that the trapped molecule is 
released at each step of the MC. The steady state probability 
distributions for this simple MC are given by 
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Since the diffusion channel is low-pass, the variations in the 
channel are slow relative to the mixing time of the MC, which 
is the time that takes for the MC to reach steady state. In 
other words, the transition time to reach the steady state is 
quite small compared to the time intervals that information 
is transmitted via the diffusion channel. Therefore, we may 
assume that the MC reaches the steady state before the 
concentration level at the receiver switches from high to low 
or vice versa. 

Let IIi be the estimator for m and denote by II i the random 
variable representing m. Again, we have: fti = YliLi Xi> 
We use the method in the previous section to obtain the op- 
timal distribution for II i to maximize the mutual information 
described in below: 



C = max /(IIi;IIi)= max [if(IIi) 

ir^Oi) /n^Tri) 
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However, the difference is that m is limited to [0, j^] as p 
varies in the range [0,1]. Consequently, II i cannot take an 
Arcsine distribution as in the ideal receptor case. 

One approach for finding the distribution of II i is to 
use Lagrange multiplier to maximize © with the constraint 
ji/(i+q) y^^) — i Since the Lagrange method is in- 
tractable, instead, the capacity achieving distribution for II i 
is numerically calculated by Blahut-Arimoto algorithm [fill 
for several values of q. The result is shown in Fig. |6(a)[ 
We observe that the distribution of II i is not symmetric 
and its range would decrease by increasing the parameter q. 
From the distribution of IIi, we can determine the optimal 
distribution for P; which in turn will determine the required 
distribution of the transmitter output. The distribution for P 
is shown in Fig. |6(b)| We observe that the distribution of P 
does not depend considerably on the value of q. Since the 
transformation from p to tt\ is one-to-one and deterministic, 
we can calculate the capacity of the ligand receptor receiver 
by computing the mutual information between II i and II i. 
The numerical results for different number of receptors and q 
are shown in Fig. |U We observe that the capacity is increased 
when q in decreased. It can be explained by the fact that, 
except for the case of q = which MC would be meaningless, 
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Fig. 6. The capacity achieving distribution (a) tti and (b) p for different values of q. 



the range of 7Ti increases when (7 decreases and hence, the gives an upper bound for overall information transmission rate 
mutual information is increased. of it. 



IV. Conclusion and Discussion 

This paper is motivated by the fundamental limits on 
the transmission rate of information in the diffusion-based 
molecular channel with ligand receptors. The problem may 
be divided into two components. The first one corresponds 
to the fundamental limit on the information transmission rate 
due to the inherent memory in the diffusion channel which 
was discussed in [ 8 ] . The second bottle-neck is at the receiver 
due to limits on concentration estimation by the receiver. 

In our previous work, we had studied a discrete diffusion 
channel model in which the transmitter maps the binary bits 
to low and high concentration levels. Furthermore, due to 
the memory of the channel, the transmitter was assumed to 
change the duration for transmitting each bit depending on the 
immediate sent bit. This dependency results in four symbols 
with different durations and the capacity of this channel was 
calculated by using Shannon's model for discrete noiseless 
system. In this paper, we presented two models accounting for 
the uncertainty in estimation of concentration of molecules by 
the ligand receptors. Using these models, we obtained the rate 
limit due to the ligand receptors. The capacity achieving distri- 
bution for input was obtained for our model which showed that 
in order to achieve higher rates, the transmitter should mostly 
produce concentrations which are either low or high. It means 
that medium levels of concentrations should be avoided which 
supports our assumption for binary transmission. 

Therefore, the results from the analysis of the receiver are 
consistent with the transmitter model we used in [8| to com- 
pute the capacity of the diffusion channel. The combination 
of results for these two components, i.e., the diffusion channel 
and the ligand receptor, gives a full model for diffusion-based 
molecular communication and cascading the two components 
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